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With the use of the symmetry-unrestricted cranked SHF method in the 3D 
coordinate-mesh representation, a systematic search for the SD and HD rotational 
bands in the N=Z nuclei from ^^S to **Cr has been done, and SD and HD solu- 
tions have been found in 223, se^r, ^Oq^^ 44^^;^ -^^ 36 ^r, 40^^^^ 44rpj^ 48^;.^ 

respectively. The SD band in ^''Ca is found to be extremely soft against both 
^ ■ the ajcially symmetric (yso) and asymmetric (yai) octupole deformations. Possible 

Ch ' presense of SD states in neutron-rich sulfur isotopes from *^S to ^^S has also been 

investigated, and deformation properties of neutron skins both in the ground and 
SD states are discussed. 



1. Introduction 

Quite recently, superdeformed(SD) rotational bands have been discovered 
in ■^^Ar, ""^Ca, and "^^Ti.^'^'^ One of the important new features of them 
is that they are built on excited 0+ states and observed up to high spin, 
in contrast to the SD bands in heavier mass regions where low-spin por- 
tions of them are unknown in almost all cases. In this talk, we shall first 
report results of the symmetry-unrestricted, cranked Skyrme-Hartree-Fock 
(SHF) calculations for these SD bands. The calculation has been carried 
out with the use of the fully three-dimensional (3D) , Cartesian coordinate- 
mesh representation without imposing any symmetry restriction. ^'^'^ The 
computational algorithm is basically the same as in the standard one '' 
except that the symmetry restrictions are removed. For comparison sake, 
we also carry out the standard symmetry-restricted calculations imposing 
reflection symmetries about the {x,y)-, {y,z)- and (z, a;)-planes. By com- 
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paring these results, we can clearly identify effects of reflection symmetry 
breaking in the mean field. 

We shall next present results of the cranked SHF calculation for SD 
bands in the neutron-rich sulfur isotopes near the neutron drip line. These 
nuclei are expected to constitute a new "SD doubly closed" region asso- 
ciated with the SD magic numbers, Z = 16 for protons and A'' ~ 30- 
32 for neutrons. An interesting theoretical subject for the SD bands in 
these neutron-rich region is to understand deformation properties of neu- 
tron skins. We shall discuss on this point. 

2. ^"Ca region 

Figure 1 shows deformation energy curves evaluated by means of the con- 
strained HF procedure. Solid lines with and without filled circles in these 
figures represent results of unrestricted and restricted calculations, respec- 
tively. In both cases, we otain local minima corresponding to the SD states 
for 32s, 36 Ar, ^Oc^ and ^^Ti in the region 0.4 < j32 < 0.8. (The local 
minimum in *^Ti is triaxial so that it is not clearly seen the 7 = sec- 
tion.) The local minima in ^^S and ^^Ar involve four particles in the fp 
shell, while those in ""^Ca and ^^Ti involve eight particles. In addition to 
these SD minima, we also obtain local minima in the region /32 > 0.8 for 
*°Ca, ^^Ti and ^^Cr. These minima involve additional four particles in the 
single-particle levels that reduce to the gg/2 levels in the spherical limit. 
Somewhat loosely we call these local minima ''hyperdeformed (HD)." The 
HD solution in *°Ca corresponds to the 12p-12h configuration. 

We notice in this figure that the crossings between configurations in- 
volving different numbers of particles in the fp shell are sharp in the re- 
stricted calculations, while we always obtain smooth configuration rear- 
rangements in the unrestricted calculations. The reason for this different 
behavior between the unrestricted and restricted calculations is rather easy 
to understand: When the parity symmetry is imposed, there is no way, 
within the mean-field approximation, to mix configurations having differ- 
ent number of particles in the fp shell. In contrast, as illustrated in Fig. 2, 
smooth crossover between these different configurations is possible via mix- 
ing between positive- and negative-parity single-particle levels, when such 
a symmetry restriction is removed. Octupole deformation parameters /^s 
are plotted as functions of in the lower portion of Fig. 3 for the case 
of '"'Ca. We see that values of are zero near the local minima, but rise 
in the crossing region. This means that the configuration rearrangements 
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Figure 1. Deformation energy curves as functions of the quadrupolc deformation P2 
calculated at / = by means of the constrained SHF procedure with the SkM* inter- 
action. The axial- asymmetry parameter 7 is constrained to be zero. The curves for 
different nuclei are shifted by 20 MeV in order to facilitate the comparison. Solid lines 
with and without filled circles represent the results obtained by the unrestricted and re- 
stricted versions, respectively (see the text). The notations /"g*" and indicate 
the configurations in which the /7/2 shell (/p shell) and the shell are respectively 
occupied by n and m nucleons. 

in fact take place through paths in the deformation space that break the 
reflection symmetry. 

Excitation energies of the SD and HD bands in ''"Ca calculated by 
using difl'erent versions (SIII, SkM*, SLy4) of the Skyrme interaction are 
compared with the experimental data ^ in the left-hand portion of Fig. 4. 
The SD band is slightly triaxial with 7 = 6°-9°(8°-9°) and it terminates 
at / c:i 24 for the SIII (SkM*) interaction. (In the case of '"'Ti, the shape 
is more triaxial with 7 = 18°-25° and 13°-19°, and the SD band termi- 
nates at 7 ~ 12 and 16 for the SIII and SkM* interactions, respectively.^) 
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Figure 2. Level crossing between single-particle levels having opposite parities. Config- 
uration rearrangement can take place by breaking the reflection symmetry in the mean 
field. The r^Ysi matrix element between levels with the asymptotic quantum numbers 
[321 1] and [200^] is large, and significantly contributes to the tendency toward the 
non-axial Y31 octupole deformation discussed in the text. 




Figure 3. The upper part: Deformation energy curve for *''Ca (same as in Fig. 1). The 
lower part: Octupole deformation fSs obtained by the unrestricted SHF calculation with 
SkM*, plotted as a function of /Si- 



February 8, 2008 23:8 WSPC/Trim Size: 9in x 6in for Proceedings 



proc 



5 



cranked SHF ( cranked SHF ) - ^{{AJ^f) 




5 10 15 20 25 30 5 10 15 20 25 30 35 

Angular Momentum / 



Figure 4. Comparison between tlie excitation energies of the SD and HD bands in '"'Ca, 
calculated by using different versions of the Skyrme interaction and the experimental 
data (filled circles). Solid, dashed and dashed-dotted hnes indicate the results with the 
SIII, SkM* and SLy4 interactions, respectively. Results with and without including 
the zero-point rotational energy correction are shown in the right- and left-hand sides, 
respectively. 

One may notice that the excitation energy of the SD band-head state is 
overestimated. We have evaluated the zero-point rotational energy correc- 
tion, ((AJa;)^), as a function of the rotational frequency a>rot- Here, 
A Ja; = Jx — (Jx) and J*^^-* denotes the dynamical moment of inertia defined 
by J(^) = dI/du>rot- Excitation energies including this correction are shown 
in the right-hand portion of Fig. 4. We see that the calculated excitation 
energies are significantly improved by including this correction. 

Let us examine stabilities of the SD local minimum in '^''Ca against 
octupole deformations. Figure 5 shows deformation energy curves as func- 
tions of the octupole deformation parameters /33m(|™| = 0,1,2,3) about 
the SD minimum. We immediately notice that the SD state is extremely 
soft with respect to the f3zo and /?3i deformations. 

Quite recently, Imagawa and Hashimoto have carried out a selfconsistent 
RPA calculation in the 3D Cartesian- mesh representation on the basis of the 
SHF mean field, and they have obtained, for the SHI (SkM*) interaction, 
a strongly collective octupole vibrational mode with = 1~ at about 
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Figure 5. Deformation energy curves (measured from the energy at f3:j = 0) as functions 
of the octupole deformation parameters /33m(|"i| = 0,1,2,3), calculated for '*''Ca by 
means of the constrained HF procedure with SkM*. The quadrupole deformation pa- 
rameters are fixed at (82 = 0.6 and 7 = 8°. One of the /33m{\'>n\ = 0, 1,2,3) is varied 
while the other Pzrn's are fixed to zero. 



1.1 (0.6) MeV excitation from the SD band head.^ Thus, it is extremely 
interesting to search for negative-parity rotational bands associated with 
the non-axial K"^ = 1~ octupole vibrational modes built on the SD yrast 
band. 



3. region 

Figure 6 shows deformation energy curves for neutron-rich sulfur isotopes 
from ^^S to ^^S, which indicates that the SD local minima is deepest at ^°S. 
As shown in Fig. 7, this result is common for the SHF calculations with 
the use of SIII, SkM* and SLy4 interactions. Thus, the neutron SD shell 
structure seems to be slightly modified from that known in the Zn region 
with N ^ Z, where N ~ 30-32 arc the SD magic numbers. 

Figures 8 and 9 show the neutron and proton density profiles for the 
ground and the SD states, respectively. We see that deformed neutron skins 
are present in both cases. These calculations are done with use of a small 
mesh size of 0.25 fm. In order to examine deformation properties of these 
neutron skins in more detail, we have made a least square fitting to the den- 
sity distribution along each principal axis direction with the Woods-Saxon 
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Figure 6. Deformation energy curves as functions of /32 calculated at / = by means of 
the constrained SHF procedure with the SLy4 interaction. The axial-asymmetry param- 
eter 7 is constrained to be zero. The curves for different nuclei are shifted by 10 MeV to 
fa<;ilitate the comparison. Solid lines with and without filled circles represent the results 
obtained by the unrestricted and restricted versions, respectively. The notations /"a'" 
indicate the configurations in which the /7/2 shell and the 99/2 shell are respectively 
occupied by n and m nucleons. 

function. The half-density radii and the diffuseness parameters extracted in 
this way for the SD state are Usted in Table 1 . We see that the neutron skin 
is formed mainly due to the difference in the diffuseness between protons 
and neutrons (rather than the difference in the half-density radius). It is 
interesting to note that the proton diffuseness parameter along the major 
axis is significantly smaller than that along the minor axis. 

The presence of the deformed neutron skins may be detected through 
excitation spectra of these nuclei. Thus, search for soft K"^ — 0^ and 1^ 
(dipole + octupole) vibrational modes in unstable nuclei with deformed 
neutron skins seems especially interesting. Note that octupole modes will 
be mixed with dipole modes in deformed nuclei. For studying these modes, 
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Figure 7. Comparison of deformation energy curves for ^"S obtained by using different 
versions of the Skyrme interaction. The curves for different versions are shifted by 10 
MeV. 

Table 1. Root- mean square radii ( -Rrms ), quadrupole deforma- 
tion parameters { 132 ), half-density radii for the major and minor 
axes ( R^^2 ' ^1/2 surface diffuseness parameters along the 
major and minor axes ( , ), evaluated by the SfiF method 
with SLy4 for the SD state in ^'^S. The proton and neutron con- 
tributions are separately shown together with their differences. 



SD 


-Rrms 


ft 


^1/2 


nx 
^1/2 






neutrons 


4.15 


0.59 


5.65 


3.09 


0.72 


0.61 


protons 


3.75 


0.70 


5.76 


2.77 


0.38 


0.47 


differences 


0.41 


-0.11 


-0.12 


0.32 


0.34 


0.14 



we need to develop the SHF-Bogoliubov -|- quasiparticle RPA approach such 
that we can take into account deformation, pairing, and continuum effects 

simultaneously. We can further envisage to go beyond the quasiparticle 
RPA by means of the selfconsistent collective coordinate method.^ 
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Figure 8. Neutron and proton density distributions in the ground state of ^"S, calcu- 
lated by ttie unrestricted SHE witli SLy4 and with a mesii size of 0.25 fm. Lcft-liand side: 
equi-density lines with 50% and 1% of the central density in the {x, z) and (x, y) planes 
are drawn. Solid and dashed lines indicate those for protons and neutrons, respectively. 
Right-hand side: Density distributions along the major and minor axes are drawn by 
thin-solid lines with filled circles. The least-square fits of them with the Wood-Saxon 
function are also shown by solid lines. The former densities are shifted up by 0.05 fm~^ 
to facilitate the comparison. 




Figure 9. The same as Fig. 8 but for the SD state. 
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4. Summary 

With the use of the symmetry-unrestricted cranked SHF method in the 3D 
coordinate-mesh representation, we have carried out a systematic theoret- 
ical search for the SD and HD rotational bands in the N—Z nuclei from 
to 48Cr. We have found the SD solutions in ^^S, ^SAr, ^"Ca, ^4x1, and 
the HD solutions in ^^Ar, ^°Ca, "^^Ti, ^*Cr. Particular attention has been 
paid to the recently discovered SD band in ^°Ca, and we have found that 
the SD band in ^°Ca is extremely soft against both the axially symmetric 
(Yso) and asymmetric (Yai) octupole deformations. Thus, it will be espe- 
cially interesting to search for negative-parity rotational bands associated 
with non-axial K'^ = 1~ octupole vibrations built on the SD yrast band. 

We have also discussed possible presense of SD states in sulfur isotopes 
from ^^S to ^^S, which are situated near the neutron drip line. An interest- 
ing subject in this region is the appearance of deformed neutron skins both 
in the ground and SD states. The presence of the deformed neutron skins 
may be detected through excitation spectra of these nuclei. Thus, search 
for new kinds of soft K"^ = 0~ and 1~ (dipole -|- octupole) vibrational 
modes of excitation is challenging, both theoretically and experimentally. 
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